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Abstract 

The updated CDF measurement of the forward-backward asymmetry ^4pB m the top quark produc- 
tion pp — > tt at Tevatron (with yfs = 1.96 TeV) shows a deviation of 2a from the value predicted 
by the Standard QCD Model. We present calculation of this quantity in the scenario where colored 
unparticle physics contributes to the s-channel of the process, and obtain the regions in the plane of 
the unparticle parameters A and du , which give the values of the ^4fb and of the total ti production 
cross section compatible with the present measurements. 
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Due to C-parity invariance, it is known that the forward-backward asymmetry (FBA) of 
top quark pair production at Tevatron vanishes at leading order (LO) |l| in the standard 
QCD model (SM). The inclusive non-zero charge asymmetry can be induced by (i) radiative 
corrections to quark-antiquark annihilation and (ii) interference between different amplitudes 
contributing to gluon-quark scattering qg -> tiq and qg — > tiq . Within the SM, this 
leads, at the Tevatron with yfs = 1.96 TeV, to the nonzero but relatively low prediction 



Af B (SM) = 0.050 ±0.015 



3, 



(1) 



Measurement of any significant deviation from this SM prediction could be attributed to 
the new physics effects. 

When DO Collaboration published the first measurement on the FBA in top-quark pair 
production in the pp laboratory frame with 0.9 fb _1 of data, an unexpectedly larger FBA 
value was indicated Qj. By using 1.9 fb _1 5], the CDF Collaboration observed the asym- 
metry (at y/s = 1.96 TeV) to be 



Af B = 0.17 ±0.08 



Ap B = 0.24 ±0.14, 



(2) 



in the pp frame and it frame, respectively. The updated CDF result with luminosity of 3.2 
fb _1 , in the pp (lab) frame, is 6] 



Af B (exp) = 0.193 ± 0.065(stat) ± 0.024(syst) . 



(3) 



As seen in Eq. fl3]), the large value of the FBA of top-quark is not smeared by the statistics. 
Inspired by the 2a deviation of the observed value of the top quark FBA from the SM- 
predicted value, several possible solutions have been proposed and studied by authors in 
Refs. |7-23|. 



In general, the top-pair production by the new physics could be through s-, t- and u- 
channel and the situation depends on the property of the new particle. No matter to which 
channel they contribute, the extensions of the SM in the framework of particle physics, such 
as axigluon js], [ll|, Z' (9], W 10], diquarks 13], etc., have some drawbacks. For instance, in 
order to explain the observed top-quark FBA value, one has to introduce unimaginably large 
flavor-changing couplings in the t- and u-channels. The couplings in the s-channel could be 
as large as the strong gauge coupling of the SM. However, beside the serious constraint from 
the invisible production of a new resonance, the sign of the top-quark couplings has to be 



chosen opposite to that of the light quarks in order to get the correct sign of the FBA value. 
In order to avoid the aforementioned problems, we study in this work the top quark FBA in 
the framework of unparticle physics which is dictated by the scale of conformal invariance. 

An exact scale-invariant "stuff" cannot have a definite mass unless it is zero. Therefore, 
in order to distinguish it from the conventional particles, Georgi named the "stuff" unpar- 
ticle 



24. 



25] . It was found that the unparticle has a noninteger scaling dimension du and 



behaves as an invisible particle 



24J (see also 26]). Further implications of the unparticle 



to collider and low energy physics are discussed in Refs. [27H29]. We will adopt three as- 
pects of unparticle physics in order to present a possible explanation of the aforementioned 
large value of the FBA. Firstly, if we take the protecting symmetry to be exact, then, due 
to the unique character of indefinite mass, no visible resonant unparticle will be produced 
in the pp collisions. Secondly, by utilizing the noninteger scale dimension, the differential 
cross section for ti production could be enhanced without fine-tuning the large couplings of 
unparticle and quarks. Finally, to match the interaction structure of the SM, the considered 
scale invariant stuff (unparticle) is a vector boson and carries color charges 30] ; it has chiral 
couplings to quarks and its representation in SU(3) C belongs to color-octet. 

Since there is no well established approach to give a full theory for unparticle interactions, 
we study instead the topic from the phenomeno logical viewpoint. In order to escape the 
large couplings from flavor changing neutral currents (FCNCs), the couplings of unparticle 
to quarks are chosen to be flavor conserving. Hence, we write the interactions of colored 
unparticle with quarks as 

C = q(g q v l» + 9 q Al^)T a qO a J , (4) 

where g x = Z^/A^ -1 and x = V an d A. Here, Z| is the dimensionless coupling and the 
index q denotes the quark flavor, Au is the scale at which the unparticle is formed, and 
{T a } = X a /2 are the SU(3) C generators (where A a are the Gell-Mann matrices) normalized 
by tr(T a T b ) = 5 ab /2. The power dy — 1 is determined from the effective Lagrangian of 
Eq. fl4]) in four- dimensional spacetime when the dimension of the colored unparticle { 



is taken as du- By following the scheme shown in Ref. 3l], the propagator of the colored 



'//< 
u 
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vector unparticle is written as 



d A xe~ ik - x 

-%Cy 



(0\TO;(x)O b u (0)\0) 

6 ab r o 2(d, 



(-P 2 



16 



\3-d u 



u 



d, 



u 



(5) 



with 



C 



A 



v 



du 



2 sin duTr 

16tt 5 / 2 T{d u + 1/2) 



(6) 



(2tt) 2 ^ T(d u - l)T(2d u ) ' 

After introducing the interactions of unparticle with quarks and the virtual unparticle 
propagator, we can now calculate the ti pair production at the quark level. Using Eqs. (OS) 
and ©, the scattering amplitude for qq — > ti by unparticle exchange in the s-channel is 



,3-du 



P 2 9vv 



-PuPu 



(7) 



where flavor q denotes the light u and d quark, and p = p q +Pq = Pt+pi- The t-channel does 
not contribute, due to flavor- conserving vertices Eq. (j3J). The equations of motion imply 
qj>q = and qfry^q = —2m q q^ 5 q. Thus, the factor 2(d u — 2)/{du — 1) in the propagator 
is associated with the light quark mass and is negligible. Consequently, the scattering 
amplitude combined with the SM contributions is given by 

A = Asm + A u 



9l 



q-i^T a qt^T a t 



+ 



sC 



v 



i3—du 



-iir(3-d u ) - ( q 



Q {qW + 9aW) T a qi (g v r + <?a7 m 7 5 ) T a t 



(8) 



with g s being the strong coupling of the QCD SM and s = (p g + Pq) 2 = (pt + Pt) 2 - For 
explicitly showing the differential cross section in ti invariant mass frame, we choose the 
relevant coordinates of particle momenta as 



1,0,0,±1 



Pt,t = — (l,±&sin0,O, ±&cosi 



(9) 
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with 01 = l — Amf/s. The polar angle 9 is the relative angle between outgoing top-quark and 
the incoming g-quark. The spin and color averaged amplitude- square is straightforwardly 
obtained as 

1 1 



\A[ 



A 



2 2 iV 2 
(N 2 g - 1) 
16iV 2 



4(47ra s 



(l + /3 2 cos 2 £ + 4^T) 



+ 8(CV47ra s ) cos7r(3 — d u ) 



<7«<&(l + #cos :4 + 4 



2g q A g t A f3 t cos6 



+ As' 



sC 



(9v) 2 ((9 q v) 2 + (ti) 2 ) l + A 2 cos 2 <? + 



4m? 



+ (9a) 2 ({9 q vf + isfif) l + A 2 cos 2 ^- 



Ami 



^9v9v9a9a^ cos 9 



(10) 



As a consequence, the differential cross section for qq — > tt process as a function of 9 in tt 
frame is found to be 



dcos9 



N2 ° 9 1 A ( (4vra s ) 2 (l + /3 2 cos 2 (9 + 4^? 



2CV(47ra s ) cos7r(3 — d, 



u 



2^ A cos 9 



+ 



g 2 cy 

S3— dh. 



(9v) 2 ((g q v ) 2 + (g q A ) 2 ) i + A 2 cos 2 ^ 



4m? 



+ (G#) 2 + Gfl) 2 ) i + A 2 cos 2 ^ 



Am? 

* + 8^<^Acos0 



In the s-channel, only the terms linear in cos# will contribute directly to the forward- 
backward asymmetry From Eq. ( flTT) . the relevant effects are associated with g\g\ and 
gygyg\g\i in which the former is from the interference between unparticle and SM while 
the latter is from the contribution of unparticle itself. In both terms, we see clearly that the 
axial-vector couplings are the essential to generate the FBA. 

To obtain the hadronic cross section from the parton level, we have to consider the 
convolution with the parton distribution functions. Thus, the differential cross section at 
the hadronic level is 

da(pp — > tt) j 1 

Jx2=0 J xi=0 



dcos9 



dx 1 dx 2 fi{x 1 )f j {x 2 ) q 

o cos 9 



(12) 



where fa (fj) is the parton distribution function of the parton (qj) in the proton (antipro- 
ton), the angle 9 represents the angle between the three-momentum of the produced t quark 



and the three-momentum of the proton p of the quark q) in the lab system (center of 
mass system of pp). The sum over (i,j) in Eq. f fl2|) is over all parton pair combinations 
qq = qiCjj for the scattering process q^j — > £t = it, d, s). 

In the following, all the unprimed kinematic quantities are in the lab system, and all the 
"hatted" kinematic quantities are in the center of mass system (CMS) of qq (<& CMS of tt). 
Taking into account the relations p q = X\p v and p q = X2Pp in the lab system, considering 
the four-momentum conservation in the scattering qq — > tt in the qq CMS, and relating the 
lab and qq CMS quantities via the corresponding boost relations, the following relation can 
be obtained between the angle 9 and its qq CMS analog 9 = 9(9, X\, x 2 ): 



cos (9(9,x 1 ,x 2 j) = -=-77 ■ ^ — 1 2Q , ^r<! - {x\ - x\) sin 2 9 

J p t l(x 1 + x 2 y - cos 2 9(xi - x 2 y\ 



+4 cos 9 



2 

x\x\$l -(xx - x 2 ) 2 sin 2 9 



1/2 



s 

where (3 t is the aforementioned quantity involving s = (p q + p q ) 2 = Xix 2 s 



(13) 



m 2 



Pt = Pt(?xX2) = \l -■ (14) 

XiX 2 S 



The relevant independent kinematic quantities in the integration ffl2|) are all lab-related: 
X\, x 2 , and 9. On the other hand, Eq. ( fl~3|) shows that the qq CMS-related angle 9 is a 
function of the aforementioned three independent quantities x±, x 2 , and 9. The quantity 
da q ^ tt / d cos 9 appearing as integrand in Eq. (112j) is obtained directly from Eqs. ( TTTT) and 
( |T3l) by applying the derivatives at fixed x\ and 2:2 



da^ tl (9, x u x 2 ) da^ tl dcos {Q( d > x ^ x ^ 



(15) 



dcos 9 d cos 9 dcos 9 

where the partial derivatives d/dcos9 are at fixed X\ and x 2 . 

The total hadronic cross section a(pp — > tt) and the corresponding forward-backward 
asymmetry are then obtained by the corresponding integrations of the expression (1121) in 
the lab frame 

a(pp^tt) = r 1 dco S 9 da( f^ , (16) 
J_i dcos 9 

A% = { f dcosO *«* -> *> - f dcosO M f ^ ) /a(pp -> tt) . (17) 
\J dcos 9 J _ x dcos 9 J 



6 



Another physical observable of experimental interest is the invariant mass distribution 



da/dM a 



32] where M t | = (p t + Vt) 2 = X1X2S 



da{pp - j. tt) = 2 Afo 



M t 2 f /s x l 



2_^M x i)fj{ x 2) / dcosO- 



1,3 



dcosO 



The integrations in Eqs. (116j) - (118l) are performed across the kinematically allowed regions, 
i.e., such that (3 t and cos# are real. 

After having presented formulas for the three physical observables, we can now numeri- 
cally investigate the unparticle contributions to top-quark pair production. At first, in order 
to reduce the number of free parameters, we assume that the colored vector unparticle is 
flavor blind, i.e. g v — g\ — g v — g q A — g. Then, the remaining unknown parameters appear- 
ing in the physical quantities are: g = X/A^~ l , and the scale dimension du. This means 
that in such a case we have only two independent parameters A and du, both dimensionless 
quantities, and we can fix the scale An formally to an arbitrary value. We will set it equal 
to An = 1 TeV. We will see that, unlike the situation in the axigluon model (si, ll| in which 



-g\ is necessary to get the positive sign in FBA, in unparticle physics the flavor-blind 
and chirality-independent couplings are enough to fit the data. 

Further, it is necessary to consider the measurements of at least two of the aforementioned 
three observables, namely o~(pp — > tt) and A^ B , Eqs. (|T6|) - f|T7|) . in order to restrict the area 
of the parameters A and du. The value of the ti production cross section a(pp —> tt) was 
measured by the CDF Collaboration 33 1 



a(pp ti) cxp = 7.50 ± 0.31 (stat) ± 0.34 (syst) ± 0.15 (th) pb 

= 7.50 ± 0.48 pb . (19) 



On the other hand, the SM prediction is a(pp — > tt) SM = 6.73~toJil pb 34J, which includes 
the contributions from the tree-level, the next-to-leading order in a s , and the next-to-leading 
in threshold logarithms (LO+NLO+NLL). In the specific case of using the CTEQ6.6 parton 
distribution functions 35] (which we use), the central value for the SM prediction goes 
slightly down to a(pp — > tt) SM = 6.61 pb, Ref. 34] (Cacciari et ai, 2008) when the top 
quark (pole) mass is taken to be m t = 175 GeV. 

The other measurement is the aforementioned FBA value Eq. ([3]). The NLO effects in the 
QCD SM give nonzero FBA value 0.050 ± 0.015, Eq. (TT]). However, in our calculations, the 




FIG. 1: ti production cross section (the lower) and top-quark FBA (the upper figure) as a func- 
tion of the scale dimension du, where the solid, dashed, dotted and dash-dotted lines represents 
A=1.4, 1.6, 1.8, 2.0, respectively. The band in the plot represents the measured values with la 
uncertainties. 

SM amplitude is the tree-level amplitude [rescaled accordingly in order to obtain a(pp — > 
ti) SM = 6.6 pb, see below], which gives A^ B = 0. Therefore, we will regard the A^ B as 
calculated according to Eq. ffTTl) [using Eqs. f|T5|) . ffl~2]) and ffTT]) ] to be responsible for the 
deviation of the experimental from the SM FBA value 

Ap B = Af B (exp) - Apg(SM) = 0.143 ± 0.071 , (20) 

where the uncertainties (±0.065, ±0.024, ±0.015) were added in quadrature. 

In our calculations we use the CTEQ6.6 parton distribution functions, the value m t = 175 
GeV for the t quark (pole) mass, and for the QCD coupling the value a s ~ a s (m t ) rs 0.11. 
With such values, we obtain the tree-level <r(SM; tree) ~ 4.85 pb. We use for the SM 



amplitude the rescaling factor V±36 in order to obtain cr(SM) = 6.6 pb, which is according 
to Ref. 34| (the second entry: Cacciari et a/., 2008) the central value of cr(SM) when 
m t = 175 GeV and CTEQ v6.6 is used for the parton distribution functions. 

The physically interesting regime for unparticle physics is 1 < du < 2, and A ~ 10°. 
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FIG. 2: As Fig. [IJ but as a function of the parameter A. The solid, dashed, dotted and dash-dotted 
lines represents du=l.l, 1.15, 1.2, 1.25, respectively. 

We calculated Afb, cr and da/dMu, scanning over the parameter regions < A < 3.5 
and 1 < du < 2. The numerical results are presented in Figs. [T]l5j We show the tt 
production cross section and [A P p B — t4^ b (SM)] as a function of d u (A) in Figs. [Q (J2J), 
where the solid, dashed, dotted and dash-dotted lines represents A = 1.4,1.6,1.8,2.0 
{du = 1.1,1.15,1.2,1.25), respectively. Figure [3] is the central result of our calculations. 
It shows the region in the dy-\ parameter plane which simultaneously fulfills the exper- 
imental constraints (|T9|) and (|20|) . This region lies between the two dotted and simul- 
taneously between the two dashed lines. The central measured values o pa 7.5 pb and 
A%(exp) ~ ^fb( SM ) ~ °- 14 are achieved at A = 2.05 and d u = 1.28. In Fig. El we scanned 
over the free parameter space in finite steps AA = 0.1 and Ady = 0.01. 

In Fig. H] we present the average values of da jdM t i in eight different Mfj-intervals ("bins") 
as used by the CDF measurement 32|. The presented results are for: (a) the QCD SM case 
(solid line: A = 0; a{tt) = 7.5 pb); (b) the "central" case (dash-dotted line; A = 2.05 
and du = 1.28, giving a{tt) = 7.5 pb and FBA of Eq. (1201 equal 0.14); (c) another, more 
"marginal" case (dashed line; A = 1.70 and du = 1.175, giving a(tt) = 7.01 pb and FBA of 
Eq. (1201) equal 0.178). The CDF measurements are the points with vertical lines. All our 
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FIG. 3: Contours for a(pp — > tt) (dotted) and A^ B — A^ B (SM) (dashed) as a function of dy 
and A. The numbers in the plot denote the lower and upper bounds of each observable with la 
uncertainties of the data: a = 7.50 ± 0.48 pb; ^4pg as in Eq. ([201) , 



results (including the QCD SM) are above the CDF measurements, with the exception of 
the first two bins M t i < 450 GeV. The deviations could plausibly be ascribed to two main 
uncertainties j36|: 

• (i) The chosen scales of renormalization (hr) and factorization (hf) for which the usual 
possible values could be taken between mt/2 and 2m t . Here we adopted = /if = m t . 

• (ii) The M^-dependent NLO effects which include the NLO parton distribution fuction 
(PDF). Here for simplicity we just use a M t i~ independent scale factor value of K=1.36 
(i.e., the factor \fK = y/1.36 for the tree-level SM amplitude Asm) to fit the tt SM 
production cross section with LO calculations. 



For a detailed analysis of the various uncertainties see Ref. [36] . Nonetheless, most of our 
results for da jdM t i are at least marginally compatible with the CDF results, within 2a. 

In order to make a more detailed inspection, in Fig. we present the M t f-restricted 
forward-backward asymmetries, i.e., those calculated by the expression f lTT|l where the phase 
space in the numerator and the denominator is restricted by M t i < M^ ge (the quantity 
Apg W ) or by M t i > M^f ge (the quantity A F g lgh ). These asymmetries were calculated for 
the two aforementioned choices of parameter values (A, dy), and are compared in the Figure 
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FIG. 4: da/dM t i as a function of invariant mass of top-pair M t ^, where the solid, dash-dotted 
and dashed lines represent the SM result and colored unparticle with (A, du) = (2.05, 1.28) and 
(1.70, 1.175), respectively. The vertical bars are the data from CDF measurement with an inte- 



grated luminosity of 2.5 fb , Ref. 



with the CDF measured values [37J subtracted by the (LO+NLO+NLL) SM values [38]. 
This subtraction is needed for comparison with our results, for the same reason as in the 
(unrestricted) A l FB of Eq. (1201) . 



4b = ^FB(exp) - Ap B (SM) (X = low, high) . (21) 



We see that the experimental uncertainties are very large, especially for Apg lgh at 
600 GeV or higher. Nonetheless, the central experimental values appear to suggest the fall 
of A^ h function of M^ ge at M^ gc > 600 GeV. It is interesting that a model involving 
axigluon exchange does give such a fall for at least one (benchmark) choice of parameters 
(with: g q A = —g^) [ill. On the other hand, our model does not show such a behavior. 
This issue remains inconclusive because of: (i) the aforementioned very large experimental 
uncertainties of ^4pJ h at high M^ ge ; (ii) the severely restricted phase space at high M^ ge . 
Namely, our simplified approach of rescaling the tree-level SM (QCD) amplitude by a fixed 
factor ( \[K = Vl-36) for all M t i values becomes increasingly unreliable when M^f ge increases 
in A F g lgh , because the phase space becomes so severely restricted. 

In conclusion, we investigated whether colored flavor-conserving unparticle physics can 
explain the measured forward-backward asymmetry value for the ti production at the Teva- 
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FIG. 5: Restricted forward-backward asymetries Apg W and ^4pg lgh as functions of the threshold 
("edge") M a values, for (A, d u ) = (2.05, 1.28) (circles) and (1.70, 1.175) (squares). Included are 



also the corresponding CDF measured values 
values 38j, as bars with triangles. 



371 ] (their 8th and 9th figure) subtracted by the SM 



tron; the latter measured value shows 2<r deviation from the QCD SM value. With a natural 
assumption of quark flavor-blind and chirality-independent interactions to unparticle, our 
calculations indicate that the aforementioned unparticle contributions can explain this devi- 
ation. We found an area of the (two-)parameter space of the unparticle physics which gives 
the results compatible with the measurements of the forward-backward asymmetry and of 
the total cross section for the ti production at the Tevatron. The resulting values of the 
differential da/dM t i cross section and the M t i- restricted forward-backward asymmetries are 
only marginally compatible with the measured values. 
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